Dark matter from dark energy-baryonic matter couplings 
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We present a scenario in which a scalar field dark energy is coupled to the trace of the energy 
momentum tensor of the baryonic matter fields. In the slow-roll regime, this interaction could 
give rise to the cosmological features of dark matter. We work out the cosmological background 
solutions and fit the parameters of the model using the Union 2 supernovae data set. Then, we 
develop cosmological perturbations up to linear order, and we find that the perturbed variables 
have an acceptable behavior, in particular the density contrast of baryonic matter grows similar to 
that in the ACDM model for a suitable choice of the strength parameter of the coupling. 



PACS numbers: 98.80.-k, 98.80.Cq, 95.35.+d 



I. INTRODUCTION 

Currently, the most successful model of cosmology 
that we have is ACDM. It is constructed in order to 
match a wide variety of modern cosmological observa- 
tions that have stunned the physicist community in the 
last decades. Among them, there are precision mea- 
surements of anisotropies in the cosmological microwave 
background radiation Q , baryon acoustic oscillation 
and Type la supernovae [5|-l7|. All these observa- 
tions point out that our Universe at present is dominated 
by a cosmological constant -dark energy- and there are 
about 5 times of some unknown nonbaryonic, dark mat- 
ter, over the baryonic matter which is well understood 
by the standard model of particles. 

Nonetheless, some objections to the ACDM model ex- 
ist, both theoretical @, Q and observational (see for ex- 
ample Thus, alternative proposals have appeared 
in the literature giving rise to the idea that dark energy 
varies with time. Scalar fields have attracted special at- 
tention, mainly quintessence [lll - [l3| . among some other 
alternatives Q. 

Recently there has been a lot of interest in studying 
couplings in the dark sector species (T^ - [23| . This is in 
part motivated by the fact that until today we can only 
extract information of these components through gravi- 
tational interaction, a feature that has been dubbed Dark 
Degeneracy [24- 28] . Specifically, we can define the energy 
content of the dark sector, and in fact the dark sector it- 
self, using the Einstein field equations 
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(1) 



where G^,y comes from the observed geometry of the Uni- 
verse and T°^'^ from its observed energy content. In this 
sense the dark sector refiects our lack of knowledge. The 
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easiest way -mathematically and conceptually- to accom- 
modate these ideas into the observed history of the Uni- 
verse is decomposing T^"''*'' in two species, dark energy 
and dark matter. 
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But this decomposition is not unique. In fact, the dark 
sector could be composed by a large zoo of particles and 
complicated interactions between them. Or, it could be 
even just one unknown field. To accomplish this last 
possibility we note that in the ACDM model the equation 
of state parameter of the total dark sector is given by 
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where i index the dark matter component and the cos- 
mological constant, and in the last equality we have used 
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In order to mimic this model 
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with just one dark field, we must have Wdark — Weff- 
Any fluid with equation of state parameter equal to Wef / 
will produce the same expansion history of the Universe. 

On the other hand, some string theory-inspired models 
of dark energy share the peculiarity that scalar fields, like 
the dilaton, couple directly to matter with gravitational 
strength. To have cosmological infiuence at present, these 
fields must be nearly massless, leading to long-range fifth- 
forces and to large violations of the equivalence principle. 
Some mechanisms have been proposed in order to avoid 
this unacceptable behavior, as the Damour-Poliakov ef- 
fect [2^, [so], in which the interaction is dynamically 
driven to zero by the expansion of the Universe, and as 
the chameleon mechanism [3l|, [s^ , where the mass of the 
scalar field has an environment density dependence, be- 
coming very huge in overall high density regions, such as 
those in which Einstein principle equivalence and fifth- 
force search experiments are performed. 

In this work we follow the both lines of thought out- 
lined above. To this end we consider an interaction La- 
grangian between a nearly massless scalar and the trace 
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of the energy momentum tensor of the ordinary matter 
fields given by 



'9A[<f>)T. 



(4) 



This type of couphng has been investigated in a cos- 
mological context in [33 35]. This interaction has some 
attractive properties, the field does not couple to the elec- 
tromagnetic field and in this sense is dark; also, it does 
not couple to relativistic matter and then do not affects 
the success of the early Universe cosmology, although it 
could couple to the inflaton field. As we shall see, the 
coupling could also give mass to the field. 

At a more fundamental level we can consider, in a first 
approximation, a fermionic matter free fields with energy 
momentum tensor Tjji, = —iil>l iiBvip and trace given by 
T^^^ — —iTp^ff^dfiip = —m^^ip, where in the last equality 
we have used the Dirac equation. The coupling then be- 
comes a Yukawa- like one and gives mass to the fermions. 
If we chose correctly the function A{ip) we can also inter- 
pret this result as that the interaction has given mass to 
the field cj). 

This paper is organized as follows: in Section II, we 
present the details of the general theory; in Section III, 
we derive the background cosmology equations and we 
choose a specific model that mimics the ACDM model, 
then we numerically obtain the cosmological solutions; in 
Section IV, we work the theory of linear perturbations; 
finally, in Section V, we present our conclusions. 



II. THE GENERAL THEORY 



The action we consider is 




167rG 



(5) 



where the coupling of the scalar field with the trace of the 
energy momentum tensor of the ordinary matter fields is 
given by (|H) . There is a difficulty here because the matter 
fields not only appears in the matter action, but also in 
the interaction action. Thus, we have to define the trace 
of the energy momentum tensor of baryonic matter as 
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(6) 



To solve this redundant definition we will be specific 
and we will work with a perfect fiuid of dust in the matter 
sector, the Lagrangian is Cm — —P\/—g, where p is the 
energy density of the fiuid in its rest frame [l^ . The trace 
of the energy momentum tensor becomes (see Appendix, 
or for an alternative derivation see l35i) 
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Then, the total action can be written as 



(7) 
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where we have defined the function e"'-'^'' through 
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These theories reduce to those with an interaction term 
f{(f>)Cm, proposed in [2^, in the case in which the matter 
fields are perfect fluids, although the interaction vanishes 
for relativistic fluids (see Appendix). The fleld equations 
derived from action |5] are 



G^. = 8^G(T^^)-|-e"(«r^™)), (10) 
where the energy momentum tensors of the flelds are 

T^f = 0,/.-^,. - ^5A..0'>.a - g^.Vi<P), (11) 



and 
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For the scalar fleld the evolution equation is 
□ 0- V^'(0) =a'(0)e"('^V, 



(12) 



(13) 



where prime means derivative with respect to (p. This 
equation shows that, as in the chameleon models, the 
evolution of the scalar fleld is governed by the effec- 
tive potential Veff{(j)) = V{(t)) + e^^'^V: and that the 
mass of slow oscillations about the minimum is given by 
— V^/j f{(j)min), where 4>min is the value of the scalar 
fleld that minimizes the effective potential. By using the 
Bianchi identities the conservation equation for the mat- 
ter fields becomes 



-a'(0)p(g^^ + u^u^)d'^(f>. 



(14) 



From this equation we obtain the continuity and 
geodesic equations 



and 



V''ipu^)^0, 
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The geodesic equation ([161) shows that only the chan- 
ges of the scalar field perpendicular to the four velocity 
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of a particle affect its motion. In particular, if the vec- where 
tor field u'^ defines a global time, only the spatial gradi- 



ent of the field influences the particle's motions. This is 
the case of the homogeneous and isotropic cosmology, in 
which there is a family of preferred comoving observers 
with the expansion. If the field is also homogeneous and 
isotropic, these observers are affected by the field only 
through gravity. 

In the nonrelativistic limit with static fields, the 
geodesic equation for a point particle reduces to 



-V$7v-a'V0-a'y(y •V)0, (17) 



where X is the spatial position of the particle, V its three 
velocity and $7v = —5^00 the Newtonian potential. The 
last term can be dropped out because it is of the order 
(w/c)^, but we prefer to leave it there because it is a 
dissipative term that arises due to the interaction with 
the scalar field and it shows explicitly how the energy 
transfer occurs. Equation p7|) also shows that the force 
exerted by the scalar field upon a test particle of mass m 
is given by 



F4, = -mV(a(</))). 



(18) 
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If the field is in slow-roll regime we neglect the 0^ con- 
tribution and the dark equation of state parameter be- 
comes 



^dark 
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(25) 



Then, the following cosmological scenario is plausible 
without dark matter. During the radiation era, the trace 
of the energy momentum tensor of the baryonic matter 
is equal to zero and the background cosmology is as the 
standard big bang. Then, in the early times of matter era 
the energy density p is high, and thus the interaction term 
dominates in the last expression, giving rise to Wdark ~ 
0. As the time goes on, p is redshifted as ~ a~^, the 
potential becomes important and eventually dominates 
over the interaction term, and in this limit Wdark — — 1- 

To obtain an acceptable late time cosmology we have 
to impose a set of constrictions in the free functions of 
the theory. First, to obtain equation ([3]) we impose the 
constriction 



III. BACKGROUND COSMOLOGY 

Considering a spatially flat FRW metric, we write the 
background evolution equations as 



H' = 



^{\^' + V{cb) + {e'^^^^-l)p + p), (19) 

?^ + 3i7(/) + y'((?;.) + a'((/))e"('^V = 0, (20) 
p + SHp = 0. (21) 



The last equation is a consequence of (ITil) (see also the 
discussion after equation (1161) ). and it can be integrated 
to give p = poa~^, where po is the value of the energy 
density of baryonic matter today and we have normalized 
the value of the scale factor today to ao = 1. This set of 
equations is a general feature of interactions of the type 
f i(t))C,n = -/(0)/5 (see for instance [H M, H^)- The 
conservation equation suggests to interpret the interac- 
tion as a pure dark sector feature. Then we define 
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(22) 



where = 0^/2 + V{4)). The total energy density is 
then pt — Pdark + P, and the Friedman equation, = 
STrG{pdark + p) ■ The conservation equation for pdark is 



Pdark + 3i/(l + Wdark)Pdark = 0, 



(23) 
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Second, to account for the correct amount of dark mat- 
ter, here given by the interaction term, we impose the 
constriction 



- 1 



(0) 
DM 
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(27) 



Ci and C2 are general functions of the scalar field that 
we expect to have a very little variation in the slow-roll 
regime. Also, to obtain the observed accelerated expan- 
sion of the Universe, we must have 



SnGViM = A, 



(28) 



where 0o is the value of the scalar field today and A is 
the measured cosmological constant. Note that this last 
equation is not independent from the two constrictions. 

The interaction over fermions impose other constraints 
because, as we show in the Appendix, its effect is to shift 



Fol- 



the fermion masses from m to m{(f)) 
lowing [38] the abundance of light elements constrains 
the variation of protons' and neutrons' masses to be 
at most of about 10% from nucleosynthesis until to- 
day. In our model this is translated into Am(^)/m(0) — 
(e"!-^") - e"('^o))/e"(*°) < 0.1, where is the value of 
the scalar field at nucleosynthesis and (f>Q is the value of 
the scalar field today. According to this the scalar field 
has to be in slow-roll at least since that epoch. 
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Fixing the free functions 



Numerical Analysis 



In the chameleon model the potential is of the runaway 
form and the coupling with matter is an exponential one. 
The effective equation of state parameter becomes of the 
form Wdark = " (1 + 0(exp(/30) — l)a~'^)~"'^ and the argu- 
ments shown above are applicable as long as the field is 
in the slow-roll regime. Nevertheless we will follow a dif- 
ferent approach and with the only aim to mimic as far as 
possible the ACDM model we choose the free functions 
of the theory as 



With the selection of the free functions given in last 
section, the cosmological background evolution equations 
can be written as 



H = —(2'^ +2"^</>'^ +0^0— + —)' (36) 



and 



(37) 



and 



(29) 



(30) 



where e is a constant with inverse squared mass units. 
Note that with this selection the minimum of the effec- 
tive potential is always at (/) = 0, and the evolution field 
equation for the scalar field becomes 



To do the numerical analysis, it is convenient to define 
the dimensionless strength parameter of the interaction 



(3 = eMl 



(38) 



Using i/o = IQAKm/s/Mpc, A/3H^ = = 0.73 and 
equation (I35p . the required value of the scalar field mass 
as a function of /3 and C2(0o) is 



= 3.6 X 10 



C2(0o) 



Mpc"^ 



(39) 



(31) 



and it is explicit that the field has an effective environ- 
ment dependent mass 



(32) 



and that the force mediated by the field has a range 
l/nieff. Although this is a feature on chameleon fields, 
the evolution equation here is linear and then we do not 
expect to obtain the thin-shell suppression [3l|, , also 
a feature in chameleon fields. 

The two constriction equations (pS)) and (P?]) become 



= 0.315 



^2 n 



(33) 



and 
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DE 
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The equation that relates the scalar field mass with A 
becomes 



\ 4>o 



(35) 



where the reduced Planck mass is given by Mp 
(8^G)-i/2. 



We evolve the equations ([55t and (|57)) for different val- 
ues of /3, using the numerical constrictions (1551) and p4p . 




a 



FIG. 1: Shows the evohition of the equation of state parame- 
ter of the dark sector, Wdark- The dashed line is the result in 
the ACDM model. The thick oscillating curve corresponds to 
P = 1; the short-dashed, to P — 0.2; and the thick nonoscil- 
lating to /3 = 0.04. 

Figure [T] shows the results for Wdark ■ The equation of 
state that is oscillating comes from a field that is oscil- 
lating about the minimum aX (p — Q and then it is not 
in the slow-roll regime. The nonoscillating curves come 
from fields that are in the slow-roll regime. Note that 
the curve with (3 = 0.04, which corresponds to e « G, 
comes from a field in slow-roll regime. 

We plot the scale factor a{t) as a function of time. 
The results are shown in figure [H We also show the 
results for the ACDM model, and the initial conditions 
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/3 = 0.01 





0.00 0.05 0.10 0.15 



FIG. 2; Evolution of the scale factor as a function of cosmic 
time. The dashed lines are the results for the ACDM model. 
The vertical lines denote the present time. Note that the 
temporal scale is different for each case. 
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FIG. 3: Evolution of C2{4>) as a function of the scale factor 
using: /3 = 0.01, thick line; and /3 — 0.04, dashed line. 



are such that at some early time both models have the 
same amount of dark matter, here given by the interac- 
tion. Note that the scale factor grows more slowly than 
in the ACDM model. This is because the emulated dark 
matter component does not redshift with a~^, but rather 
with (jp'a~^ . The age of the Universe is, in the case of 
/3 = 0.01, about 3% larger than in the ACDM model, and 
about 10% for the case (3 = 0.04. 

Figure |3] is a plot of the function C2 {(f>) as a function 
of the scale factor for the strength parameters /? — 0.01 
and /? = 0.04, with initial value C2{(f>i) — 5. Note that in 
both cases, after an initial transient period, the functions 
become less steep. 

Figure |3] shows the evolution of the different density 
parameters as a function of the number of e-folds, N = 



FIG. 4: Evolution of the different density parameters as a 
function of TV = log a, we use /3 = 0.01. The upper line is 
^^dark; the dotted line is Qt; the dashed line, Qint', and the 
thick line, fi^. The vertical line denotes present time. 



logo. These are given by 

8nG / 1 ■ 2 1 2 , 2 1 , 2 Po \ 



Note that fidark — ^int + and fidark + i^b = 1- 



Fit to supernovae data 

Figure [5] shows the best fit to the Union 2 data set 0, 
a recent sample of 557 supernovae. We use P = 0.01 and 
the minimum value of the function is given by Ci = 
0.305 and C2{4>o) = 4.66, where 0o is the cosmological 
value of the scalar field today. The module distance /x 
and the luminosity distance are related by 

Figure [6] shows the contour plots at Icr and 2<j for this 
same fit. 

In the background cosmological solutions the model 
looks very similar as the ACDM one as far as /? ^ 0.04. 
In the next section we consider the linear perturbation 
theory of the cosmological solutions. 



IV. LINEAR PERTURBATIONS 

In the last section we have shown how the background 
dynamics can be made very similar to the ACDM model. 
This is mainly because the background fluids are affected 
by the scalar field only gravitationally, see equation ([2T|) . 
This comes ultimately from the fact that the right-hand 
side of the conservation equation projects out tem- 
poral variations of 0, with respect to the observer u'^. 
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FIG. 5: Supernovae Union 2 data and the prediction for the 
model (soUd Une). We have used /3 = 0.01, Ci = 0.305 and 
C2{(po) = 4.66. 



Union2 



^ 0.30 




FIG. 6: Contour plots at la and 2a. The central point is Ci = 
0.305 and C2(0o) = 4.66, which corresponds to a minimum 
value Xd.o.f. = 0.98. 



In the background cosmology only temporal variations 
of the scalar field are considered, and so it does not exert 
any net force over the baryons. This is not true in the 
inhomogeneous cosmology where spatial gradients of the 
fields have to be taken into account. Accordingly one 
expects that the evolution of cosmological perturbations 
of baryonic fields could differ from the obtained in the 
ACDM model. In this section we show that the model 
has an acceptable growth of baryonic matter density per- 
turbations. 

We consider scalar perturbations in the longitudinal 
gauge. The metric is given by 



-(1 + 2'^)dt^ + a^it){l - 2$)%da;Ma;^. (41) 




-0.007 L 



FIG. 7: Evolution of the initial perturbations for a comoving 
wavenumber k = 0.05Mpc~^, using jS = 0.01. For subhori- 
zon scales k ^ H, the relation S{k2)/S{ki) = ki/ki for the 
density contrast holds approximately. The dashed lines are 
the results in the ACDM model and the vertical lines denote 
present time. The background dynamics are the same as in 
figure [2l 
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The fields perturbations are given by 



to 
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FIG. 8: Tiie same as figure [7] but witli /3 — 0.04, which corre- 
sponds to e ~ G. 



Pb{t,x) 
(j>{t,x) 



pit){l + dit,x)), 



(42) 
(43) 
(44) 



where p{t) and Uq are the background energy density 
and velocity of the baryonic matter fluid, and 4>o{t) is the 
background value of the scalar field. As usual we define 
the variable 9 of the baryonic fluid as the divergence of 
the peculiar velocity of the fluid in momentum space (see 
[37j for details). 



(45) 



where ki is the i component of the comoving wave num- 
ber vector A: of a perturbation, related to the physical 
wavelength by A = I'Kajk. 

In this analysis we neglect anisotropic stresses and then 
the scalar potentials are equal, = $. The perturbation 
equations in Fourier space are given by 



+3i/(* iJ*) = 47rG - <^aH - rnl(i>Ql 



(46) 



2i.2 . 



(47) 



47rG[ -(^2* + ^oH ~ ml^bdcj)], (48) 
and the hydrodynamical perturbation equations 



S - + -9 ^ 0, 
a 



(49) 




and 



= -2TO|0o*-epo0o('5 + 2*)-h4?!)o*. (51) 

These six equations are not all independent, we use 
(|47|. (149)) . ((50)) and ((STj) to numerically evolve a set of 
initial perturbations. The background dynamics are the 
same as in figure [51 In fact, for the case /3 = 0.01, they 
yield the best fit to supernovae data in the last section. 
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The results of the numerical analysis are shown in figures 
[7] and [5] The dashed lines are the results in the ACDM 
model. 

An interesting parameter is the squared of the effective 
sound speed of the perturbations of the scalar field, = 

(52) 

For an uncoupled scalar field dark energy component, 
as quintessence, the sound speed is exactly one in the 
slow-roll approximation, preventing scalar field structure 
to form. In our model, the scalar field acts both as dark 
energy and, by the coupling to baryons, as dark matter. 
As a dark energy component, scalar field perturbations 
tend to be erased, but on the other hand the interaction 
drives to zero the speed of sound squared, due to the last 
term in the denominator of equation (j52p . allowing the 
scalar field perturbations to grow. As a result of this 
interplay scalar field perturbations are not completeley 
damped, instead initial scalar field inhomogeneities could 
be amplified or at least frozen. This behavior is shown in 
figures [7] and IH where the scalar field perturbations do 
not decay to zero as the time goes on, but they tend to 
a negative value about which they oscillate. 

The scalar field perturbations become a source of the 
gravitational potential (see equation that acts upon 
the baryonic matter through equation This is a 

feature of dark matter. However, note that the sign of 
the average of the scalar field perturbations is negative 
and, from the Poisson equation (46), they act as a repul- 
sive gravitational source. Nevertheless, the gravitational 
potential, read from figures [7] and [U is similar to the 
one in the ACDM model. This is because the baryonic 
matter energy density appears multiplied by the factor 
(1 -|- £0^/2) in equation (|46p and it contributes as an 
attractive gravitational source. On the other hand, the 
scale factor grows at slower rate than in the ACDM model 
and this yields an increase in the gravitational perturba- 
tions through the second term of the left-hand side, of 
equation P5|) . As a net result of these effects the gravi- 
tational potential is slightly stronger (^I^ more negative), 
and therefore the density contrast grows faster, than in 
the ACDM model. These can be seen in figures [7] and [H 

As it was pointed out in (39j coupled dark energy mod- 
els could suffer from unwanted severe instabilities char- 
acterized by a negative speed of sound squared. This is 
not the case in our model in which the numerical analysis 
shows that the scalar field perturbations do not grow but 
instead oscillate around a nonzero value. A more detailed 
analysis [2l| has shown that these instabilities are sup- 
pressed when the scalar field is in the slow-roll regime, 
as in our case. 



V. CONCLUSIONS 

It is possible that future observations and experiments 
reveal that there exist some interactions between the dark 
sector and the known fields of the standard model; in- 
deed, this is the hope for terrestrial experiments in de- 
tecting dark matter. The dark degeneracy allows us to 
study a wide variety of models that mimic the ACDM at 
today's observation accuracy but with a richer dynamics. 

In this paper we have presented a plausible mecha- 
nism in which an interaction between a quintessence-like 
field and baryonic matter gives rise to effects similar to 
those of dark matter. The interaction is given through 
the trace of the energy momentum tensor of the matter 
fields, and in the case that we have considered, dust mat- 
ter, it is the same as an interaction of the type f{(j))C 
The differences are manifest if we instead use a radiation 
field, e.g. photons, in this case the interaction vanishes 
-neglecting the trace anomaly (40| - and in this sense is 
dark. Because of this coupling the radiation era is the 
same as in the standard big bang and the differences 
appear once matter dominates. We have described the 
background cosmology and by fitting the parameters of 
the model we can mimic as far as we want the late time 
ACDM background solutions, from decoupling until to- 
day. Specifically we have shown the differences for the 
cases of /3 = 0.4 (e w G) and /3 = 0.01. For this last 
value of /3, we made fits to the Union 2 supernovae data 
set and found that the best fit is obtained by choosing the 
other two parameters of the model equal to Ci — 0.305 
and C2(0o) = 4.66. 

We have worked out the first order perturbation the- 
ory and shown that these models are capable to form 
linear structure without dark matter. Initial perturba- 
tions of a free scalar field dark energy are erased because 
its sound speed is equal to one, and by the fact that 
dark energy accelerates the Universe, it tends to freeze 
any matter perturbation. One of the effects of the inter- 
action, when the field is slow-rolling, is to decrease the 
sound speed. In our model this induces the scalar field 
perturbations to oscillate about a nonzero negative value 
and yield a repulsive gravitational force over the baryonic 
matter perturbations. Despite this effect, we showed that 
the baryonic matter density contrast could grow as fast 
as in the ACDM model. This is because the interaction 
enhances the gravitational strength of the baryons. This 
latter effect increases as fi does, but at the same time the 
increase of /? can break the slow-roll. Thus, suitable /3 
are found for /? < 0.04. 



VI. APPENDIX 

In this Appendix, we show how the interaction affects 
different types of matter, namely a fermion field and a 
perfect fluid. 

For a fermion Dirac field we focus on the microscopic 
interaction and we consider a fixed Minkowski spacetime. 
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The Lagrangian is £ = C^+dnt+C^ = {-■ip{ir]'^''jf_td^- 
m)'ip + A{4))T)^—ri + £0 , and the trace is given by equa- 
tion ([6]), which now reads 



where we have used the identity 5p = i(p 



T = — 



2 5{S^, + Sint) 
1= c ' 



(53) 



9ti}^)S9^'' (see for example [36i |l 
obtain the differential equation 



dT 



Taking the trace, we 



For a free fermion field, one obtains T^^^ = — iV'T/^^^V' 
—mipil) , where we have used the Dirac equation. Now,_for 

the total Lagrangian we make the ansatz T = me'^^'^^ip , with general solution 
and inserting it into equation ()53p . we obtain that 



3^(0)(1 + w)p— + (1 - 'iA{(t,))T + (1 - 5w)p = 0, (57) 
dp 



(1 — . So, the Lagrangian now reads 

C = -^A(«7^5,, - e"("^)m)V' + C4,. 



(54) 



The effect of the interaction is to shift the mass of the 
fermions from m to e"'"^^m. 

Now, let us consider in the matter sector a perfect 
fluid with an equation of state parameter w constant. 
The matter Lagrangian is Cm = —p^J—g-, where p is the 
energy density of the fluid in its rest frame [s^ . We will 
assume that the trace of the energy momentum tensor 
depends only on p, T = T{p). It also depends on the 
metric, but only through p-. Then 6T = (dT / dp)5p . 
Let us consider the action 



S^nt + Srn = / d"" xV^{Aicj))T<-"'^ -p), (55) 



performing variations with respect to the metric, equa- 
tion (O gives 



-^(0)-^p(l + w)(u^-u^ -^5^1,), (56) 



T ■ 



(1 -3w) 



1 - (1 - Sw)A{(j)) 



p + Cp 



'(l-4A)/3A(l+w) 



(58) 



where C is an integration constant that we have to choose 
equal to zero if we want for a radiation fluid T = 0. Then 
we obtain the next equation for the trace 



T 



1 - 3w 



1- {l-3w)A{(t>)' 



(59) 



For dust, w ^ and this equation reduces to (O. For 
a relativistic fluid w = 1/3 and the trace vanishes. This 
result agrees with the one obtained through an iterative 
procedure in '33, . 
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